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The Dirac structure of the nucleon self-energy in asymmetric nuclear matter cannot reliably be 
deduced from the momentum dependence of the single-particle energies. It is demonstrated that 
such attempts yield an isospin dependence with even a wrong sign. Relativistic studies of finite 
nuclei have been based on such studies of asymmetric nuclear matter. The effects of these isospin 
components on the results for finite nuclei are investigated. 
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I. INTRODUCTION 

During the last years a substantial progress has been made in the microscopic description of the bulk properties 
of nuclear matter by the inclusion of relativistic features in the so-called Dirac-Brueckner-Hartree-Fock (DBHF) 
approximation. One finds that the self-energy of the nucleon in nuclear matter contains a large attractive component 
Ss which is of the order of -300 MeV and transforms like a scalar under a Lorentz transformation. This very attractive 
contribution is compensated to a large extent by a repulsive, time-like Lorentz-vector component Eg. This partial 
cancellation between Eg and Eg leads to single-particle energies and binding energies for the nucleons, which are of 
the order of -40 MeV, small compared to the nucleon rest-mass. Based on this small binding energy it has been argued 
for a long time that relativistic effects should be small in nuclear physics. 

Flowever, considering a self-energy with a relativistic structure as just outlined in a Dirac equation for a nucleon in a 
medium of nuclear matter, one finds that Dirac spinors derived from this equation exhibit a substantial enhancement 
of the small component as compared to the Dirac spinor of a free nucleon with the same momentum. This Dirac spinor 
essentially corresponds to one for a nucleon with an effective mass m* which is the sum of the bare mass m plus the 
scalar part of the self-energy Eg. Assuming a value for Eg of -300 MeV, which is quite typical for nuclear matter at 
saturation density, it is obvious that this reduction of the effective mass can result in non-negligible effects. The matrix 
elements of the nucleon-nucleon (NN) interaction for two nucleons in nuclear matter should be evaluated employing 
these Dirac spinors modified by the nuclear medium, which means that the NN interaction is density-dependent. It is 
this density-dependence which leads to the saturation of nuclear matter in simple relativistic mean-field calculations 
within the Walecka model 

Replacing the phenomenological approach for the NN interaction used in the Walecka model by a realistic meson- 
exchange potential, one has to perform a nuclear structure calculation which goes beyond the mean-field or Hartree- 
Fock approximation and account for the effects of NN correlations. This can be done in a relativistic extension of 
the Brueckner-Hartree-Fock approximation. For that purpose one has to evaluate the G-matrix by solving the Bethe- 
Goldstone equation. For relativistic models of the NN interaction like the One-Boson-Exchange (OBE) potentials 
this Bethe-Goldstone equation corresponds to a three-dimensional reduction of the Bethe-Salpeter equation, which 
accounts for Pauli- and dispersive corrections due to the surrounding nuclear medium Q. 

The G-matrix derived from this Bethe-Goldstone equation can be analysed and decomposed into five Lorentz 
invariant amplitudes. From these amplitudes one can calculate the nucleon self-energy within the Brueckner-Hartree- 
Fock approximation and determine its relativistic structure, which means its decomposition into a scalar term Eg, a 
time-like vector term Eg and a space-like vector term E^,. This self-energy is included in a Dirac equation to determine 
the Dirac spinors for the nucleons. A self-consistent solution of these DBHF equations requires that the resulting 
spinors are used to evaluate the matrix elements of the OBE interaction and determine the G-matrix. 

Such self-consistent DBHF calculations for nuclear matter have been performed by various groups employing differ¬ 
ent models for the NN interaction H-P. All these investigations show that the relativistic effects modify the saturation 
properties of nuclear matter derived from realistic NN interaction within the BHF scheme. The density-dependence 
of the nucleon Dirac spinors yields some repulsion which increases significantly with the nuclear density. Due to this 
mechanism the Dirac effects of the DBHF approach supply a fine-tuning in the calculated energy of nuclear matter as 
a function of density, which moves the prediction for the saturation point off the well-known Coester band 0 , which 
is obtained in non-relativistic many-body calculations of nuclear matter, towards the empirical result. Brockmann 
and Machleidt actually succeeded in finding a version of the Bonn potential, which fits the NN scattering data and 
reproduces the empirical saturation point of symmetric nuclear matter using the DBHF approach |^. 
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Various groups analyzed the relativistic structure of the G-matrix and derived the three components of the self¬ 
energy, Sg, So and S^,, in symmetric nuclear matter. They find the dependence of these three components on the 
momentum of the nucleon is rather weak and that the effects of the space-like vector component St, is rather small 
as compared to the other two iiil- Therefore it seems to be justified that one takes advantage of this feature and 
derives the decomposition of the self-energy into Sg and Sq from the momentum-dependence of the single-particle 
energy e(fc) (see also below). In this way one can avoid the analysis of the relativistic structure of G, which simplifies 
the self-consistent procedure of DBHF significantly. 

After the DBHF scheme had been applied with great success to symmetric nuclear matter, it was an obvious 
extension to use the same scheme for pure neutron matter and assymetric nuclear matter with various fractions of 
proton to neutron density as well ||l^-Q. In section 2, below, we would like to demonstrate that it is very dangerous to 
apply the simplified version of the DBHF self-consistency, i.e. deducing the Dirac components from the single-particle 
energies, to asymmetric nuclear systems. We will see that this scheme, which has been used e.g. by IMl even tends 
to predict the “wrong sign” for the isovector dependence of Sg and Sq. 

Attempts have also been made to apply the DBHF approach to finite nuclei as well null- One possibility in 
this direction is the so-called Relativistic Mean Field approach with Density dependent coupling constants (RMFD) 
@0- The first step of this approach considers the Relativistic Mean Field or Dirac-Hartree approach for nuclear 
matter at a given density, adjusting an effective coupling constant for the scalar meson a and the vector meson to 
to reproduce the results of microscopic DBHF calculations at this density. In a second step these density-dependent 
coupling constants are employed in a Dirac-Hartree calculation to evaluate the ground-state properties of finite nuclei 
in a kind of local-density approximation. 

Recently this RMFD approach has been extended to asymmetric nuclear systems. Using the DBHF results of 
Engvik et al. |^, Shen et al. |]0| determined density-dependent coupling constants for the isoscalar mesons a and 
u) and the isovector mesons S and p to fit the relativistic components of the self-energy for protons and neutrons at 
various densities and asymmetries. Since, however, this analysis is based on DBHF calculations which exhibits the 
“wrong sign” in the isospin structure, one obtains coupling constants for the isovector mesons in particular, which 
are not appropriate. In section 3 we would like to correct this sign and show the effect of a proper treatment of the 
isovector mesons on the structure of finite nuclei using the RMFD approach. 


II. ISOVECTOR MESONS AND ASYMMETRIC NUCLEAR MATTER 

Our investigation of the role of isovector mesons in the effective NN interaction in a nuclear medium is based on 
the analysis of Boersma and Malfliet . They parametrize the relativistic structure of the G-matrix in terms of five 
Lorentz invariants 


G = ^T“F(“)F(“) 


( 1 ) 


with 


^(“) = 


1 , 757^, 7'", 75 


2m* 


( 2 ) 


The amplitudes T" are parametrized in terms of Yukawa functions, depending on the Mandelstam variable t as 






( 3 ) 


with effective meson masses and coupling constants gan- At a given density of nuclear matter, these coupliiw 
constants are adjusted in such a way that the antisymmetrized matrix elements of the parametrization (|^) - (|^ 
reproduce the corresponding matrixelements of the G-matrix, calculated in the rest frame of nuclear matter pS. 
More details of this parametrization can be found in . All the results displayed in this manuscript are obtained 
using the parameter set HI of Results derived from parameter sets I and H are rather similar. 

Using this parametrization one can easily determine the self-energy S*(fc) for protons and neutrons, labeled by 
index i, with momentum k 


EXfc) = IS* (fc) + 7 • kSUfc) - 7°S* (fc), 


( 4 ) 
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by solving the Dirac-Hartree-Fock equation for the NN interaction defined by the parametrization of d) - d) in a 
selfconsistent way. Inserting this self-energy into the Dirac equation for a nucleon in the nuclear medium, we obtain 


(l + S;(A))7-k+(M + El(fc)) 




Ui{k) 




( 5 ) 


Now it is convenient to eliminate and rewrite this Dirac equation into a form which only contains a scalar and a 
time-like vector component 


where 


[7 • k -f (M + lll{k)) - E* (A)70] u,{k) = e,(fc)7V(fc) 



- me; 

1 + Ej, 

1 + S* 


( 6 ) 
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We will assume that the parametrization of G, which has been determined by Boersma and Malfliet for symmetric 
nuclear matter at various densities p may also be used for asymmetric nuclear matter at the same density. With this 
assumption we ignore the fact that the G-matrix will not only depend on the density of nuclear matter but also on 
the asymmetry parameter 


a = 


P 


p 


pn + pP 


( 8 ) 


where p^ and p" denote the density of protons and neutrons, respectively. This implies that symmetric nuclear matter 
corresponds to a = 0.5 and pure neutron matter to a = 0. The assumption that the G-matrix does not depend on 
this asymmetry parameter a may not be sufficient for a very sophisticated study of asymmetric systems. It should 
be sufficient, however, for the more general remarks, which we wish to make in the present study. 

As a typical example we will now first consider the case of asymmetric nuclear matter which is defined by a baryon 
density p = 0.185 nucleon fm“^, which corresponds roughly to the saturation density of nuclear matter and an 
asymmetry parameter a = 0.35. Results for the self-energy of protons and neutrons are displayed in Fig. ^ as a 
function of the momentum k. One observes a momentum dependence of the scalar and vector components, E* and 
Eg defined in (^, which is rather weak as compared to the total value of these components. The variation of these 
quantities as a function of momentum is around 20 MeV for momenta below the Fermi momentum, which corresponds 
to about 5 percent of the total value. Therefore one may be tempted to ignore this momentum-dependence for a 
moment and interpret the mean value for these components 


Up- ^ 
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( 9 ) 


in terms of a mean field or Hartree-Fock model. The index /3 in this equation represents the scalar {/3 = s) or vector 
component (/? = 0) and kpi stands for the Fermi momentum of protons (i = p) and neutrons {i = n), respectively. 
Assuming a meson exchange model for the NN interaction, which considers the exchange of a scalar meson (cr) and a 
vector meson (w) plus the exchange of the corresponding isovector mesons S and p, the components of the self-energy 
for protons and neutrons are easily evaluated within the Hartree approximation as 


= -^ ip^s + Ps) + ^ ip^s - P^) 

Trig 

u^ = -^ip^ + pn-^ip^-pn 

mi mj, 

K = -^ip^+pn + ^{p^-pn 

mf, mi 


( 10 ) 
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The parameter and in these equations refer to the meson-nucleon coupling constants and the masses of the 
varies mesons, the proton and neutron densities are represented by and p” and the scalar densities are defined by 


Stt 



( 11 ) 


with the effective mass 


M* =M + Ul (12) 

and 

E* = yjM* + . 

Inserting these Hartree-results for the self-energy into the Dirac eq.(||) one obtains the single-particle energy in the 
Hartree approximation 


ef (fc) = ^M* + _ t/*. (13) 

For nuclear matter with p” > pP the Hartree approximation of (0) predicts more negative values for the self-energy 
components for neutrons, Uq and D", as compared to the corresponding results for the protons. Looking again at 
the momentum-dependent results for these self-energy components displayed in Fig. |I] we observe that this feature 
is reproduced by the momentum-dependent self-energy components derived from the Groningen parametrization of 
DBHF. From this figure we can furthermore observe that the isospin dependence is larger for the vector component 
of the self-energy Sg than for the scalar component. Translated into the Hartree analysis of eq.(j^) this means that 
the effective p exchange contribution to the G matrix seems to be stronger than the effects of the isovector scalar 
meson 5. This is definetely true for a realistic meson-exchange potential like the Bonn potential |^, and obviously 
remains valid for the G matrix. The difference in the isospin dependence of Eg and Eg is also responsible for the 
single-particle energies 


e.(fc) = v/(M + E|(fc))VP-E*(fc) (14) 

and leads to the result that the single-particle energies for the protons are more attractive than those for the neutrons 
(see part on the right-hand side of Fig. p. This reflects the fact that the proton - neutron interaction is more attractive 
than the neutron - neutron interaction. 

One of the main features of the dependence of the self-energy for protons and neutrons in asymmetric nuclear matter 
is displayed in the left part of Fig. ||. There we show the dependence of the average effective mass M* calculated 
according eq. (H) using the momentum averaged Ul defined in (^ as a function of the asymmetry parameter a. The 
effective mass is of course identical for protons and neutrons in the case of symmetric nuclear matter {a =0.5). The 
effective mass for the protons increases with decreasing a i.e. with an increasing fraction of neutrons. The effective 
mass for the neutrons is smaller than the corresponding M* for all values 0 < a < 0.5. 

This is just opposite to the behaviour observed by Engvik et al. (see also Fig.2 of 0) as well as Lee et 

al. 0. These investigations report effective masses for the proton to be smaller than those for the neutron. One 
could argue that this different behaviour of E* or M* is due to the different interactions: the results of have been 
derived from the Bonn potential A, whereas we consider parametrization of the Groningen potential. However, we 
would like to demonstrate that the isospin dependence of the effective mass M* reported in is wrong and that 
this mistake is caused by the fact that Engvik et al. deduce the effective mass from the momentum dependence of the 
single-particle energy, a method which has been used by other groups before 

In order to demonstrate this, we consider the momentum dependent self-energies E%(fc) deduced from the Gronin¬ 
gen potential and calculate the single-particle energy ei{k) calculated according to (|4^. In a next step we analyze 
these functions ei(fc) and deduce from these functions effective masses M* and vector potentials Uq. This can be done 
by taking at two values oi k = ki, k 2 (our choice is ki = O-lkp, ^2 = kp) and by adjusting the parameters M* and 
Uq in such a way that ( |l^ reproduces the values ei(fci) and e(fc 2 ). In this way we simulate the procedure of Engvik 
et al. and Lee et al. |13f| for the Groningen potential. 

The results of this procedure is displayed in the right part of Fig. ||. We find that the effective masses deduced 
from the momentum dependence of the single-particle energies show the same behaviour as reported by Engvik et al. 
I p^ and Lee et al. [^. This demonstrates very clearly that the isospin dependence dependence of the effective mass 
is very sensitive to the method used to determine it. Although the momentum dependence of the components of the 
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self-energy Eg ^ is rather weak on the scale of the absolute values of these quantities (see discussion of Fig. || above), 
it is obviously too strong to neglect it in analyzing the momentum-dependence of ei(fc). Therefore such an analysis 
leads to wrong results, in particular when we are interested in the isospin dependence. 

This observation led us to compare the two schemes of analyzing the Dirac structure of the nucleon also in the case 
of symmetric nuclear matter. Results for the effective mass as a function of the Fermi momentum kp of symmetric 
nuclear matter are displayed in Fig. The results are rather similar for densities around the empirical saturation 
density {kp w 1.4 fm“^). This suggests that the evaluation of the saturation properties of nuclear matter, as performed 
by Brockmann and Machleidt may not be influenced too much by improving the method to determine the Dirac 
structure of the nucleon spinors. For other densities, however, the effective mass deduced from the momentum 
dependence of the single-particle energy (dashed line) deviates significantly from the effective mass derived directly 
from Eg. In particular at smaller densities, these are the densities which are relevant for the description of finite 
nuclei with medium mass, the effective mass derived from the single-particle energies is considerably smaller than the 
corresponding value deduced from Eg. 


III. RELATIVISTIC MEAN FIELD CALCULATIONS FOR FINITE NUCLEI 


Relativistic calculations for finite nuclei have often been based on DBHF studies of nuclear matter using a kind of 
local density approximation One possibility is to consider asymmetric nuclear matter at a given density p 

and asymmetry a and determine the scalar and vector components Ug and Uq of the nucleon self-energy for protons 
and neutrons. Using the set of equations (^) one can determine effective coupling constants g^, g^, and g^. This 

means that a Dirac-Hartree calculation of nuclear matter allowing for the exchange of cr, w, 5, and p mesons with 
these coupling constants would reproduce the results of the DBHF calculation at this density and asymmetry. These 
effective coupling constants depend on density and asymmetry. 

As an example we present in Fig. ^ the density-dependence of these coupling constants, derived from the mean 
values calculated according (||). These coupling constants are represented by the solid lines and have been evaluated 
at an asymmetry defined by a =0.35. The dependence of these coupling constants on a turns out to be weak. One 
can see that these effective coupling constants exhibit a moderate density-dependence. The coupling constants for 
the isovector mesons S and p are weaker than those of the isoscalar mesons. 

Instead of using the momentum averaged self-energies one may also try to use the corresponding quantities 
Ug g derived from the momentum-dependence of the single-particle energies e(k). Using these quantities as input to 
determine the effctive meson-coupling constants, we obtain the results represented by the dashed line. It must be 
mentioned that the coupling constans g^ for the isovector mesons turn out to be negative in this analysis, which 
means that the coupling constants g themselves would be imaginary. This is a consequence of the “wrong” isospin 
dependence of these self-energy components. The effective coupling constants derived from the single-particle spectra 
also show a much stronger density dependence and the absolute values for the isoscalar mesons are close to those of 
the isovector mesons. We would like to mention again that both sets of effective coupling constants (represented by 
the solid and dashed lines in Fig. ^ are determined for the same parameterization of the Groningen G-matrix. The 
only difference is that the solid line is derived from the self-energies directly, while the dashed lines are obtained from 
the attempt to derive these self-energies from the single-particle spectrum. 

Shen et al. fl^l determined effective coupling constants in the same way as we just outlined using as input data 


the results of Engvik et al. |12|. This means that they start from components U* g which are derived from the single¬ 


particle spectrum e(k) calculated for the Bonn potential. It is remarkable that the coupling constants derived in 
show features very similar to those represented by the dashed line above. The coupling constants g^ are negative for 
the isovector mesons (note that the isovector components in eq.(6) of |l^] have a wrong sign) and all effective coupling 
constants show a strong density dependence, very similar to the dashed lines in our Fig. 

Considering a fixed asymmetry parameter a one may now consider the various coupling constants, which for the 
asymmetry assumed just depend on the total baryon density and solve the Dirac-Hartree or relativistic mean field 
equations assuming a local density approximation for the meson-nucleon coupling constant. Details of the procedure 
have been presented in j^. Note that we do not include the effect of rearrangement terms ||l7|| . 

Results for such Dirac Hartree calculations with density dependent meson nucleon coupling constants are listed 
in Tab. |. As examples we consider the ground-state properties of the nuclei and Three different models 
are considered, (a): a pure a - oj model, the density dependent coupling constants for these two mesons are derived 
from the results of symmetric nuclear matter, (b): a cr, w, g - d model in the limit of symmetric matter, i.e. a = 
0.5. The density dependent coupling constants for these mesons are obtained by analyzing the proton and neutron 
self-energies for asymmetric nuclear matter, extrapolated to the symmetric case of a = 0.5. (c): a cr, w, g - <5 model 
for asymmetric nuclear matter at a = 0.36, which corresponds the proton fraction of For each of these models 
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(a) - (c) the coupling constants derived either from the relativistic structure of the self-energy (results listed in the 
left part of the table) or from the momentum dependence of the single-particle energies (results in the right part). 

One can see from these results that the two methods to determine the density dependence of the effective meson 
coupling constants lead to quite different results even for the isospin symmetric system The difference in the 

calculated binding energy is as large as 0.6 MeV per nucleon. The deviation of the second scheme from the first one is 
even larger for the non-symmetric nucleus Note that the differences between the two methods, around 0.9 MeV 
per nucleon, is larger than the effect of the isovector mesons: the difference between the a — to model and the models 
including isovector mesons is only around 0.3 MeV per nucleon. 


IV. CONCLUSIONS 

The accuracy of the simple technique which tries to extract the Dirac structure of the nucleon self-energy in 
nuclear matter from the momentum dependence of the single-particle energy PJTo|-p^ has been investigated using 
the relativistic parametrization of the G-matrix of Boersma et al. [^. It is demonstrated that the analysis based 
on the single-particle spectrum fails if asymmetric nuclear systems are considered. The isovector components of the 
self-energy derived in this way exhibit even the wrong sign. The simplified analysis reproduces the structure of the 
self-energy around the saturation density but overestimates the scalar part of the self-energy at smaller densities. 
Therefore the improvement of this scheme may not have a very significant effect on the calculated saturation property 
of symmetric nuclear matter Q . A more apropriate analysis of the Dirac structure of the nuclear self-energy is required 
to derive density-dependent effective coupling constants to be used in Dirac-Hartree calculations of finite nuclei. 

This work has been supported by the “Graduiertenkolleg Struktur und Wechselwirkung von Hadronen und Kernen” 
(DFG GRK 132/2). 
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from S 

E/A [MeV] 

Rch [fm] 

from e(fe) 

E/A [MeV] 

Rch [fm] 

16^ 

cr, u: 

-6.94 

2.56 

-6.36 

2.60 


Exp. 

-7.98 

2.74 



22 0 

(7, U> 

-6.72 

2.56 

-5.83 

2.59 


a = 0.5 

-6.33 

2.59 

-5.47 

2.61 


a = 0.36 

-6.40 

2.58 

-5.35 

2.63 


Exp. 

-7.36 





TABLE I. Energy per nucleon (E/A) and radius of the charge distribution Rch calculated in various versions of the Dirac 
Hartree approach with density dependent coupling constants. The calculated energy has been corrected to account for the 
effects of the spurious center of mass motion. Details of the various approaches are discussed in the text. 
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k [fm-i] k [fm-i] k [fm-i] 


FIG. 1. Scalar (Ss, left part) and vector component (Eq, middle part) of the self-energy and the single-particle energy (e, 
see ( |l^ ), right part of the figure) for protons (solid line) and neutrons (dashed line) as a function of momentum k. Results 
are displayed for asymmetric nuclear matter with a baryon density p = 0.185 nucleon fm“® and an asymmetry parameter a = 


0.35. 



a a 


FIG. 2. Effective masses of nuclear matter (p = 0.185 nucleon fm“®) as a function of the asymmetry parametr a. The 
effective masses displayed in the left part of the figure have been derived from the momentum averaged self-energy of eq.(|^, 
while those in the right part of the figure are deduced from the single-particle energies e(fc) according eq.(p^ ). 










FIG. 3. Effective masses of symmetric nuclear matter as a function of the Fermi momentum kp- The results represented by 
the solid line have been derived from the momentum averaged self-energy Es of eq.(^. Those shown by the dashed line are 
deduced from the single-particle energies e(fc) according eq. (H). 



Density Density [fin-^] 


FIG. 4. Effective couplin g c onstants to be used in Dirac-Hartree calculations as a function of density. These coupling 
constants are derived from (flffl assuming meson masses of 550 MeV, 783 MeV, 983 MeV, and 769 MeV for a, to, 5, and p 
mesons, respectively. The solid lines are derived from the momentum averaged self-energies, while the dashed lines are obtained 
using effective masses derived from the single-particle spectrum e(fc). Note that the coupling constants gf are negative for the 
isovector mesons if they are derived from e{k). 
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